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We use analytic (current) density-potential maps of time-dependent (current) density functional 
theory (TD(C)DFT) to inverse engineer analytically solvable time-dependent quantum problems. 

In this approach the driving potential (the control signal) and the corresponding solution of the 
Schrodinger equation are parametrized analytically in terms of the basic TD(C)DFT observables. 

We describe the general reconstruction strategy and illustrate it with a number of explicit examples. 

First we consider the real space one-particle dynamics driven by a time-dependent electromagnetic 
held and recover, from the general TDDFT reconstruction formulas, the known exact solution for 
a driven oscillator with a time-dependent frequency. Then we use analytic maps of the lattice 
TD(C)DFT to control quantum dynamics in a discrete space. As a hrst example we construct a 
time-dependent potential which generates prescribed dynamics on a tight-binding chain. Then our 
method is applied to the dynamics of spin-1/2 driven by a time dependent magnetic held. We design 
an analytic control pulse that transfers the system from the ground to excited state and vice versa. 

This pulse generates the spin hip thus operating as a quantum NOT gate. 


I. INTRODUCTION 

Exact analytic solutions of the Schrodinger equa¬ 
tion have always been of great methodological interest 
as they underlie our intuitive understanding of quan¬ 
tum systems. Numerous solutions of the stationary 
Schrodinger equation are known since the early days of 
quantum mechanic^. However, for a long time there 
were very few examples of analytic solutions describ¬ 
ing an evolution of a quantum system driven by time- 
dependent external held, like solutions of the Landau- 
Zeneil^and RabP problems, or the solution for a driven 
harmonic oscillator^ clos^ related to a so called har¬ 
monic potential theorenPH^ . The interest in analytic 
solutions to time-dependent quantum problems was re¬ 
newed with the emergence of quantum computing. The 
necessity of designing quantum gates requires an ac¬ 
curate control of qubit dynamicJI^®^ and the state 
preparatiorP^^^. It has been recognized that by us¬ 
ing analytical pulses in quantum control problems one 
achieves a more robust evolution against errors and pulse 
parameterwhich explains the practical impor¬ 
tance of hnding new solvable quantum problems and the 
popularity of a few known pulses for the analytic control 
of two-level systems.^^^^^. 

In the last decade a number of new analytic so¬ 
lutions to the Schrodinger equation have been con¬ 
structed by inverse engineering time-dependent Hamil¬ 
tonians from given dynamics of state vectors. We note 
that most o f thes e studies focus on dynamics of two- 
level system^^^^^^^f^ and a few examples of three-level 
systems .1^ In the present work we propose an alter¬ 
native strategy of reconstructing time-dependent driv¬ 
ing potentials for analytically solvable quantum prob¬ 
lems. Our proposal employs the ideas and theorems 
of time-dependent density functional theory (TDDFT) 


and time-dependent current density functional theory 

(TDCDFT)|3te. 

Originally TDDFT/TDCDFT was developed as an 
extension of the static DFTpl for addressing time- 
dependent quantum many-body problems. The key 
statement underlying this approach is a so called map¬ 
ping theorem that establishes a one-to-one map from the 
time-dependent density/current to the external driving 
potential. The existence of this map implies that the 
knowledge of some properly chosen one particle observ¬ 
ables (collective variables, such as the density or the cur¬ 
rent) is sufficient to uniquely reconstruct the correspond¬ 
ing conjugated driving fields, and, therefore, the full wave 
function of the system. For a general many-particle sys¬ 
tem the density-potential map is practically never known 
explicitly. In come simple situations it can be constructed 
numerically. For example, recently a numerical fixed- 
point algorithm has been used to reconstruct a potential 
that produces a pr escribed time-dependent density in a 
model system.l^^® 

Remarkably, for one-particle systems the density- 
potential and the current-vector potential maps can be 
found explicitly in a closed analytic fornP^^®. In the 
present work we use the explicit TDCDFT maps to con¬ 
struct analytic control signals driving a system in such a 
way that the prescribed behavior of the basic collective 
variable, the current and/or the density, is reproduced. 
The time dependence of the control signal and the dy¬ 
namics of the wave function are then parametrized in 
terms of the physically intuitive observable. The ana¬ 
lytic TDCDFT maps are known both for a particle in 
the real continuum space and for discrete, lattice (e. g. 
tight binding) systems. This allows us to address, within 
a common scheme, control problems for the real space dy¬ 
namics and for dynamics of discrete systems with a finite 
dimensional Hilbert space, such as a motion of quantum 
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particle on tight-binding lattices, or the dynamics of a 
spin in the presence of a time-dependent magnetic field. 
To illustrate our strategy of inverse engineering we will 
recover the known exact solution for a driven harmonic 
oscillato]®^, and present nontrivial examples of analytic 
control for a particle on a finite ID chain and for a spin- 
1/2 (qubit) system. 

The structure the paper is the following. In Sec. [n| 
we present the general idea of reconstructing driving po¬ 
tentials for solvable problems using analytic TD(C)DFT 
maps. In Sec. |m] we use TDCDFT to construct solv¬ 
able problems for the real space one-particle dynamics. 
As a particular example we reconstruct the potential and 
the wave function generated by a density evolution in a 
form of a time-dependent rescaling of some initial dis¬ 
tribution supplemented with a rigid shift in space. The 
corresponding solution recove rs th e one for the driven 
harmonic oscillatoJ^. In Sec. llVl the formalism for dis¬ 
crete spaces is presented. In the first subsection we give 
an explicit example for controlling motion of a particle 
on atomic chain. In the second subsection the formal¬ 
ism is appied to a spin-1/2 control that is isomorphic to 
the control problem for a particle on a two-site lattice. 
Finally, in SecjVjwe summarize our results. 


II. CONSTRUCTION OF SOLVABLE 

PROBLEMS VIA TDDFT MAPS: THE BASIC 
IDEA 

In the standard ’’direct” statement of a quantum me¬ 
chanical problem the Schrodinger equation determines 
evolution of the wave function from the initial state 
in the presence of a given time-dependent external 
potential. Thus for a given initial state the Schrodinger 
equation generates a map V ^ 1-L from the space V of 
external potentials to the Hilbert space 1-L. This direct 
map is shown schematically on Fig. Unfortunately, 
the solution of the time-dependent Schrodinger equation , 
even for simplest two-level systems, practically always re¬ 
quires numerical calculations. Analytically solvable time- 
dependent quantum problems are exceptionally rare. 

To understand how TDDFT can help in finding analyt¬ 
ically solvable problems we analyze mapping between dif¬ 
ferent sets of object entering this approach. All TDDFT- 
type theories rely on the existence of a unique solution 
to a special “inverse” quantum problem. That is a pos¬ 
sibility to uniquely reconstruct the driving field from a 
given evolution of the conjugated observable (such as, 
the density in TDDFT or the current in TDCDFT) and 
a given initial state. In other words, if JV is the space of 
basic observables, then the existence of TDDFT implies 
that for a given initial state there exist two unique maps 
Af V and Af ^ H which relate a given trajectory in 
the space JV of observables to the corresponding trajec¬ 
tories in V and in H. The composition of these TDDFT 
maps recovers the usual direct map V ^ 1-L generated by 
the time-dependent Schrodinger equation . 



FIG. I. (a) The direct map from a trajectory in the space V of 
potentials to the trajectory in the Hilbert space generated 
by the time-dependent Schrodinger equation for a given initial 
state '00- (b) TDDFT mappings between trajectories in the 
space Af of observables (densities or currents), space V of 
potentials and the Hilbert space H. For a given initial state, 
by choosing a desired time evolution of the density/current we 
can reconstruct the driving potential and the wave function. 


In general for many-particle systems the solution of the 
inverse problem is even more difficult than the solution of 
the usual Schrodinger equation . In fact, mathematically 
construction of the TDDFT maps is equivalent to solv¬ 
ing a certain nonlinear quantum many-body problenP^. 
However there special situations when the inverse prob¬ 
lem possesses a simple analytic solution. These situa¬ 
tions cover, in particular, generic driven one-particle dy¬ 
namics both in the real space and on lattices with some 
mild restrictions on allowed initial states and the behav¬ 
ior of observables. For those cases the TD(C)DFT maps 
Af ^ V and Af ^ % can be found explicitly in the an¬ 
alytic fornP^J®!. For example, in the case of a particle 
driven by a time-dependent vector potential A(t) the 
wave function and the vector potential are the explicit 
analytic functionals of the current density j(r, t) and the 
initial state, that is A[j,'0o] and '0[j,'0o]- By construc¬ 
tion the wave function and the potential, obtained in such 
a way, are connected by the Schrodinger equation with 
the proper initial condition. Hence by assuming different 
space-time distributions of the observable we can gen¬ 
erate infinitely many solutions to the Schrodinger equa¬ 
tion , where the potential in the Hamiltonian and the so¬ 
lution are expressed analytically in term of the prescribed 
observable. In this setup the space of observables plays a 
role of the parameter space, while the TD(C)DFT maps 
provide us with the analytic parametrization formulas for 
the Hamiltonian and the solution. 

In the next sections we present the explicit framework 
for designing solvable one-particle problems in continuous 
and discrete spaces and illustrate our strategy of inverse 
engineering with several nontrivial examples. 
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III. RECONSTRUCTION OF THE REAL SPACE 
POTENTIALS 

In this section we illustrate our strategy of generat¬ 
ing the exact solutions for a simpler and more familiar 
case of a single quantum particle in the real space. Let 
us consider an electron in the three-dimensional space 
subjected to a time-dependent external electromagnetic 
field. It is convenient to use the temporal gauge in which 
the electric E and magnetic B fields are related to the 
vector potential A(r, t) as followJ^ 

B = V X A; E = -dtA. (1) 

Given the vector potential A(r, t) and the initial state 
'0o(r) = |'0o(r)|e'^°^^^ the wave function '^(r,t) is ob¬ 
tained by solving the following Schrodinger equation 


all analytically solvable and the solutions are given by 
Eqs. pbl) and (|4^. 


Now we turn to a more specific situation when the 
system is driven by a longitudinal electric field E at zero 
magnetic field, B = 0. In the absence of the magnetic 
field V X A = 0 and therefore Eq. (4a) implies that the 


curl of the velocity also vanishes, V x v = 0. The ve¬ 
locity of a one-particle system driven by a potential field 
is must be potential. In this case it is natural to use 
the Coulomb gauge (V • A = 0) and express the electric 
field as a gradient of the scalar potential U(r, t), that is 
E = — VU. By applying the standard gauge transforma¬ 
tion to Eq. ([^ we obtain the Schrodinger equation of the 
following form 


idt'ipir,t) = ( 


(6) 


= ^(-iV-A(r,i))V(r,i). (2) 

Here we are interested in the vector potentials for 
which the Schrodinger equation possesses an analytic so¬ 
lution. Therefore we follow our strategy and apply TD- 
CDET maps j i-> A and j ^ ip from the time-dependent 
current j(r, t) to the vector potential and the wave func¬ 
tion. These maps can be easily found explicitly provided 
the initial state is nodeless ipo{r) ^ 0, and the current 
j(r, t) fulfills the condition 

/ V ■j(r,t) ^ IV’oWI^ = no(r), (3) 

Jta 


which ensures that the state remains nodeless in the 
course of evolution. Then A(r, t) and the time-dependent 
wave function ip(r, t) are uniquely reconstructed from the 
given currenfp^ 


xp[tpo,j] = 


<fbPo,i] = ¥>0 + 


/'dt'fZivpgl 

Jto 


(4a) 

(4b) 




(4c) 


where the density n(r, t) and the velocity v(r,t) are de¬ 
fined as follows 


n{r,t) = -f V ■ j(r,i)+no(r), 

(5a) 

Jta 


v(r, t) = —- r. 

(5b) 


One can straightforwardly check by a direct substitution 
that these formulas indeed give a solution to Eq. §. 
These equations allow us to reconstruct the driving po¬ 
tential for a prescribed evolution of the current. 

Using Eqs. with different time dependent cur¬ 

rents we can construct infinitely many Schrodinger equa¬ 
tions with time dependent vector potentials, which are 


where '0(r, t) is now the time-dependent wave function in 
the Coulomb gauge. 

Applying the gauge transformation Eqs. we find 
the mapping from the current j(r, t) or equivalently the 
velocity and density n(r, t) to the external potential 
U(r, t) and the wave function '0(r, 


Vir.t) = - —m 

2m^/n{r, t) 


f v(r',t)-dr' 

Jo 


2^v 


C(t) 

tA(r,t) = 

0(r, t)= f mv(r', t) • dr'+ C(t), 
Jo 


(7a) 

(7b) 

(7c) 


where C{t) is a time-dependent constant. Since by con¬ 
struction the velocity is irrotational the line integrals 
in Eqs. ( [7a| ) and ( [T^ do not depend on the integration 
path. Therefore we indicate only the initial and the final 
points of the path. The value C{to) at the initial time 
is uniquely determined by the initial condition, while for 
t > to the function C(t) is arbitrary and can be chosen at 
convenience, for example to fix the value of the potential 
at infinity. The presence of a time-dependent constant in 
the density-potential mapping is in agreement with the 
Runge-Gross theorenP^. The first term in Eq. ( [7a| ) is the 
Bohm potential that can be interpreted physically as an 
adiabatic potential for which the prescribed (nodeless) 
n(r, t) is the instantaneous ground state density. The 
second and the third terms in Eq. (7a) are related to in¬ 
ertia forces. These terms compensate the inertia forces 
exerted on a particle in a local non-inertial frame moving 
with the velocity v(r,t). As a result in this co-moving 
frame the density stays stationary and equal to the initial 
density distribution. In the original frame the velocity- 
dependent contribution appears as a deformation of the 
adiabatic potential, which is aimed at supporting the pre¬ 
scribed density in the case of arbitrary fast evolution. 

There is an important difference of the present con¬ 
struction and the explicit current-vector potential map¬ 
ping of Eqs. Equations 0 state that given the 
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density n and the corresponding velocity v, the exter¬ 
nal potential V and the wave function ip can be found 
analytically. However the density and the velocity are 
not independent variables as they have to be consistent 
through the continuity equation 

t) = -V ■ [n(r, t)v(r, t)]. (8) 


corresponds to the instantaneous ground state of the 
shifted and rescaled potential a~‘^VQ {{r — yq)/ a) with 
the ground state energy a~‘^EQ. Therefore the Bohm 
potential can be represented as 


VVn(r,^) ^ _}_Y ( r-ro(0 \ 
2my^n(r, t) ^ \ a{t) ) 


a^{t) 




( 11 ) 


The complication comes from the requirement of irrota- 
tional velocity, V x v = 0 which implies the velocity field 
of the form v = V4>. Because of this condition there is 
no a simple and universal analytic relation between the 
observables entering 0- Such a relation can be found 
only for ID systems, or if we assume a ID inhomogene¬ 
ity of the observables. In higher dimensions our ability 
of constructing solvable quantum problems is limited by 
the possibility to solve analytically a classical hydrody¬ 
namics problem of reconstructing the density from the 
velocity or vice versa for an irrotational flow. Below we 
present a simple example of such a reconstruction. 


A. Exact solution generated by a time-dependent 
scaling of observables 


The final results for the potential and the wave function 
generated by the velocity 0 [or equivalently by the den¬ 
sity (10)] take the following form 


V{Y,t) = T^o f -—-mfo -r 
\ a J 

Tn a, ,9 ^ X 

-(r-ro), (12a) 

2 ol 

>*(M) = (I2b) 

, m a . 2 

i) = IT - (r - ro) +mro-r 
2 a 

~ [ + (12c) 

Jo ^ ^ 


Let the evolution starts from the ground state -00 of 
a potential Vb(r) with the ground state density no(r) 
and the energy Eq. The simplest irrotational velocity 
field v(r, t) for which Eq. 0 can be solved analytically 
is a linear function of coordinates with time-dependent 
coefficients 


v{r,t) = X{t)[r-ro{t)]+ro{t). (9) 


This veiocity corresponds to rigid motion of a fluid sup¬ 
plemented with a uniform expansion/compression rela¬ 
tively to the origin moving along the trajectory r = ro(t). 
The expansion/compression scaling factor is related to 
the parameter X{t) as, a{t) = This interpretation 

is confirmed by soiving the continuity equation ^ with 
the veiocity of Eq. 0- The corresponding soiution for 
the density takes the form 


n{r,t) 




no 


r-i'o(t) \ 

a{t) J ’ 


( 10 ) 


which indeed corresponds to a rescaied density moving 
aiong the trajectory r = ro(t). The assumed initiai con¬ 
ditions, n(r, to) = '^ 0 ( 1 *) and v(r,to) = 0, are fuifliied if 
the time dependent parameters A and tq have zero vai- 
ues and zero time derivatives at the initiai time, that is 
A(to) = 0, ro(to) = 0, A(to) = 0 and fo(to) = 0. 

Now we can insert the prescribed observabies, Eqs. 0 
and ( p^ , into Eq. 0 to reconstruct the corresponding 
potentiai and the wave function. 


To caicuiate the Bohm potentiai entering Eqs. (7a) we 


make use of the fact that no(r) is the ground state density 
of the potentiai Vb(r) with the energy Eq. This implies 
that the shifted and rescaled density n(r,t) of Eq. (10) 


Obviously, the first term in Eq. (12a) is the adiabatic 


potentiai. The other two terms describe two types of in¬ 
ertia forces - the usuai iinear acceieration force (the sec¬ 
ond term) the inertiai force reiated to a time-dependent 
deformation. 

In the speciai case of rigid motion, a = 1 or A = 0, only 
a linear acceleration inertial correction survives, so that 


the potential of Eq. (7a) simplifies as V(r,t) = Vo(r — 
Yo) — mYo • r. This potentiai rigidiy transports a quan¬ 
tum system aiong a given trajectory without any reshap¬ 
ing of the initiai density proflie. It is worth noting that 
in this particuiar case our soiution to the Schrodinger 
equation is not limited to one particle and can be triv¬ 
ially generalized to a system of any number of interacting 
identical particles. Indeed, the solution generated by a 
spatially uniform velocity held v(r,t) = Yo(t) can be ob¬ 
tained by the transformation to a uniformly accelerated 
reference fram^. Since the relative motion of particles is 
unaffected by this transformation the above potential will 
transport the center of mass while keeping unchanged the 
quantum state for the relative motion. It is absolutely 
obvious that if the initial state corresponds to that of 
the harmonic potential, our solution is identical to the 
harmonic potential theorenP^^. 

One can also easily see that the analytic solution of 
the Schrodinger equation for a harmonic oscillator with 
a time-dependent frequency uj(t) and a driving force f(t]P 

idti/>(r, t)= “ f(i) ■ V’(r, t), 

(13) 

is a particular c ase o f our Eq. (12). Assuming Vb(r) = 
^mcjQr^ in Eq. ( |12a[ ) we And that the reconstructed po¬ 
tential V (r, t) coincides (up to irrelevant constant) with 
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the potential in Eq. (13), where 



(14) 

a 

f(t) = m ^-I*g -1- mrg. 

\a^ a J 

(15) 


From these two equations we observe that the center of 
mass position rg is the solution to the Newton equation 
for a driven harmonic oscillator 


mro + muj‘^{t)rQ = f(t). (16) 

Hence in this particular case the solution of Schrodinger 
equation for the driven quantum oscillator is expressed 
in terms of the solution for the classical driven oscillator, 
which is the main observation made in Refs. m and El 


IV. INVERSE ENGINEERING OF SOLVABLE 
QUANTUM PROBLEMS ON A DISCRETE 
SPACE 


In this section we describe and illustrate our general 
reconstruction strategy for lattice systems. In this case 
we use the maps for a generalized lattice-TDCDFTS^ to 
inverse engineer analytically solvable one-particle prob¬ 
lems (or problems isomorphic to one-particle dynamics 
on a lattice]P^. 

Our starting point is the Schrodinger equation for the 
wave function '0i(to) describing a particle on an M-site 
lattice with time-dependent complex hopping parameters 

Tij, 


M 

idti’iit) = - ; Tu = 0, (17) 

i=i 


The link current Jij{t) and the on-site density n^(t) = 
|'0i(t)p are connected by the lattice continuity equation 
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Since the link current and the link kinetic energy are, 
respectively, antisymmetric and symmetric with respect 
reversing the direction of the lattice link, Jij = —Jij and 
Kij = Kij^ the combined complex observable Qij is a 
Hermitian matrix, Qij = 

Given the complex current Qij (t) and the initial state 
ipi{to) = ^ 0, the complex hopping Tij{t) 

and the wave function can be ex¬ 

pressed explicitly as functions of Qij and ipi{to) 


Tij{t) — 


Qij{t) 




\Mt)\ = . li’iitoW - [ di', 

\| j 

(t)+ C 


^i(t) = ipi 


-dt'. 


(21a) 

(21b) 

(21c) 


These formulas provide us with the analytic lattice- 
TDCDFT map from the observable to the conjugate driv¬ 
ing potential and the corresponding solution of the time- 
dependent Schrodinger equation . Using this map we can 
construct infinitely many analytically solvable problems 
generated by different time-dependent Hermitian observ¬ 
ables Qij{t). 

Below we will give two examples which illustrate the 
possibility to analytically control quantum dynamics 
time in a discrete space. 


A. Dynamics of one particle on a ID chain 


where indexes i and j take values on the lattice sites 
indicating the position in the discrete space, and Tij = 
Tj- to have a Hermitian Hamiltonian. In the Schrodinger 
equation 0 we adopted a temporal gauge in which the 
scalar on-site potential and, possibly, a magnetic field 
enter via the phase of the hopping parameters Ty(f)P. 
For generality we also allow a time-dependent hopping 
rate \Tij{t)\. 

In the generalized lattice-TDCDFT of Ref.|43]the com¬ 
plex hopping Tij{t) plays a role of a driving potential. 
The corresponding observable in this approach can be 
called a ’’complex current”!^ 

Qij{t) = 2Tij{t)i3l{t)^ljj{t). (18) 


In this subsection we use our approach to manipulate 
the on-site density of a quantum particle on a finite tight- 
binding chain. Let us consider a particle on an atomic 
chain with a nearest neighbor hopping parameters 
of of fixed amplitude \Ti^ij^i \ = Tg. The dynamics of the 
system is described by Eq. 0. Since for ID systems 
only scalar (on-site) driving potentials are allowed, one 
can always gauge transformation the Hamiltonian to the 
form with real hopping parameters = Tg and the 

real on-site potential Vi{t)^^ In the new gauge, which 
is the lattice analog of the Coulomb gauge, the time- 
dependent Schrodinger equation reads: 

idti^i = -Tg(V^i+i -h i^i-i) + Vi{t)il3i. (22) 


The real part of Qij (t) is equal to physical current Jij on 
the lattice link connecting sites i and j, while the real its 
real part Kij represents the kinetic energy on the link 

(19) 


In the following we assume for definiteness that the evo¬ 
lution starts from the ground state -^g of the chain. 

The equations for the observables (18) and (19) for 
Qij ,Jij and Kij remain the same except in the right 
hand side of Eq. (18) the hopping Tij and the density 


Qij — T^ij '^Jij' 
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matrix pij need to be replaced by their counterparts in 
the Coulomb gauge. 

The map (21) from the complex current Qij to the 
hopping Tij and the wave function in the new gauge 
is transformed to an analytic map from Qij to the on-site 
potential Vi and the wave function 


M 


s( 

j=i 




2ni 


2n 


d+i 




AT^nitij+i J ’ 

i’iit) = (23b) 

^i{t) = + f (23c) 

Jtn 


The important point is that in the considered physical 
situation with the fixed hopping amplitude the above for¬ 
mulas are not sufficient to reconstruct the potential from 
the given dynamics of observables. The fixed value of the 
hopping amplitude sets an upper bound on allowed val¬ 
ues of link currents. As a result not all possible Q become 
physically allowed, or 'T’-representable in the TDDFT ter¬ 
minology. In fact, from the definition of Eq. ( p!8| ) we find 
that the modulus of physically allowed Qij (t) is bounded 
from above 


|Qy(i)|2 = 4Tini{t)njt) < 4TI 


(24) 


Formally the condition of the fixed hopping amplitude 
reduces the dimension of the space M of observables. In 
the present case this restriction can be taken into account 
by expressing Jij (t) and Kij (t) in terms of on-site density 
ni{t). Firstly, in ID we can solve the continuity equation 
([2Q|) to get the link current 


Ji, 

i+i = 

i=i 

(25) 

Secondly, we express 
Eqs. (24) and (25) 

Kij{t) in terms of ni{t) 

using 

Ki^i^i = 

4:Tq Tliflij-l ^ ^ h jf ^ , 

1=1 

(26) 


where the sign is determined by the sign of Kij{^ at 
the initial time through the given initial state Fi¬ 


nally, by inserting Eqs. (25) and (26) into Eq. (23a) we 


obtain the explicit analytic formulas for the reconstruc¬ 
tion of the on-site lattice potential and the correspond¬ 
ing wave function from a given time-dependent density 
distribution. These form ulas correspond to the maps of 
the lattice TDDFTEIEII, Jt is interesting to note that 
exact solution proposed recently in Ref. [30] for a driven 
two-level system is, in fact, based on the above lattice 
TDDFT maps for a particular case of a two-site lattice. 


Let us now demonstrate how this map works in prac¬ 
tice by constructing a potential that a produces pre¬ 
scribed evolution of a density. Consider a particle on an 
atomic chain with II sites and a positive hopping con¬ 
stant To = I and assume that the dynamics starts from 
the ground state of the chain with zero on-site potential. 


?/)j(0) = -Ip: 




1 . /7ri. 


(27) 


We will construct the driving potential which generates 
the following two-stage evolution: (z) On the first stage 
for 0 < t < ti the system evolves from the ground state 
of Eq. ^7\ to a state with a homogeneous density distri¬ 
bution rii = I/I I; (ii) On the second stage for ti < t < t 2 
the homogeneous density distribution shrinks to the cen¬ 
ter of the chain and by t = ^2 concentrates at site 6 with 
a Gaussian envelope, rii oc . The the required 

time evolution of density ni{t) for this two-stage process 
is the following, 

( ^S{t/ti) + {1 - S{t/ti))\ipf\^ 0<t<ti 
mit) = { (28) 

i A/'(t)exp[-5( ^^)(i-6)2] ti<t<t 2 

where Af{t) is the normalization factor 


-v(() = (i; 


o ^ 12 — 11 


)(i-6)" 


-1 


(29) 


i=l 


Here S{x) is a smooth step-like function which start from 
zero at X = 0 and reaches unity at x = I. For the reason 
that will be clear later we choose a function which has a 
zero first and second derivatives at x = 0,1. Specifically 
here we use the following smooth step function which 
satisfies the above conditions conditions 


S{x) = X - sin(27rx). 

27r 


(30) 


The time dependence of on-site densities rii {t) defined by 
Eq. (28) with ti = 3 and t 2 = 12 is shown on Fig. 
Each line on the figure shows the prescribed evolution 
of the density on a particular site. At t = 0 the system 
is in the ground state (27), then it goes gradually to 


the homogeneous distribution at t = 3. Afterwards the 
density starts shrinking and finally at t = 12 it reaches a 
bell shaped Gaussian centered at the middle site. 

The analytic representation for the corresponding driv¬ 
ing potential can now be found immediately by inserting 
Eq. (28) into Eq. (23a) where the link current Jij and 


the kinetic energy Kij are given by Eqs. (25) and (26), 
respectively, and plugging in to the equation for the 
on-site potential Vi (|23a|). Since for the initial ground 


state Kin(O) is negative the minus sign must be chosen in 
Eq. (26)® Fig. shows the on-site potentials for the first 
stage of the dynamics, 0 < t < 3. Each curve represents 
the time dependence of the potential for a particular site. 
Similarly, Fig. shows the driving potential for the sec¬ 
ond stage of the time evolution. 
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FIG. 2. (Color online) The time evolution of the on-site densi¬ 
ties determined by Eq. (28). Each line represents the dynamic 


of the density in a particular site. The time duration for the 
first stage is 3 units, = 3 and for the second stage is 9 
units, t 2 = 12. Dots and their envelops indicate the density 
distribution in the initial t = 0, intermediate (t = 3), and the 
final (t = 12) states. 


EIG. 4. (Color online) On-site potentials (23a) as functions 
of time for the second stage of the evolution, 3 < t < 12. The 
evolution for 3 < t < 7, is zoomed in the magnified box. 


those times and therefore the current Jij and its first 
time deri vative also vanish. By setting Jij = Jij = 0 in 
Eq. (23a) we find the potential of the form 



EIG. 3. (Color online) On-site potentials (23a), for the first 
stage of the evolution, 0 < t < 3, as a function of time. 
The on-site potentials for all sites are zero at the beginning. 
At t = 3 all Vi except those for the boundary sites reach 1, 
while the potentials for the two ending points stay zero. Color 
online. 


Our reconstructed potential has one interesting prop¬ 
erty. On the first stage of the evolution the potential 
shown on Fig. [^drives the system from its ground state 
at t = 0 an takes it into a state with a homogeneous den¬ 
sity at t = 3 which is also the ground state of the system 
with the instantaneous potential Vi(t = 3). Similarly, 
for the second stage of the evolution the driving poten¬ 
tial shown Fig Intakes the system from ground state and 
brings it into a new ground state with a Gaussian enve¬ 
lope. At the first glance this behavior looks surprising as 
the dynamics is by far non-adiabatic. The explanation 
is, however, simple. The system at t = ti and t = ^2 is 
in its instantaneous ground state because the first and 
second derivative of the step function S (30) are zero at 


^ V V^i+1 J 

This is a lattice analog of the Bohm potential that corre¬ 
sponds to the ground state potential for a given instanta¬ 
neous density. Therefore our construction can be used to 
make a fast transfer of a system between different ground 
states. 


B. Reconstruction of a driving magnetic field for a 
spin-1/2 system 


As the last example we inverse engineer an analytically 
solvable Schrodinger equation for a driven spin-1/2 sys¬ 
tem, or, equivalently, a generic two-level system. This 
formalism can be used, for example, to control the state 
evolution of a spin-1/2 using the time-dependent mag¬ 
netic field. A similar problem has been addressed re¬ 
cently in Ref. [311 Below we apply our general reconstruc¬ 
tion strategy based on the lattice TDCDFT mapping. 

Assume a spin in the initial state |'0o) subject to a 
time dependent magnetic filed B(t). The time dependent 
Schrodinger equation for the state vector |'0(t)) reads 

idt\^{t)) = -B(i) ■ S |■^/)(^)), (32) 


where S is the spin-1/2 operator. 

By a gauge transformation one can always eliminate 
^-component of the magnetic filed and therefore reduce 
the problem to solving the Schrodinger equation with the 
magnetic field in the x^-plane 





where the ’’complex” magnetic field B is 


B={B,-iBy)/2, 


(33) 


(34) 




























and = ('0(t)| t) and = {'ip{t)\ 1) are the 

projections of the state vector on the eigenstates of Sz- 
Equation (33) is identical to the Schrodinger equa¬ 
tion 0 for a two site lattice where the spin indexes 
t and I label the sites, and B{t) is the complex hopping 
parameter. Therefore we can directly apply the lattice- 
TDCDFT maps of Eqs. ( [2T| to reconstruct the driving 
magnetic field and the wave functions 
from given dynamics of the complex observable Q{t) = 
K{t) iJ{t) defined in Eq. (18). In the present case the 
reconstruction formulas reduce to the form 


= ^l^t4(0)pT j ImQit') dt', (35a) 

mW = m(o)+r 

Jo 




m = 


Q{t) 




(35b) 

(35c) 


Equations ( [3^ provide an analytic parametrization of 
the driving field and the wave function in terms of a given 
trajectory in the two-dimensional space J\f = {K, J} of 
observables. The point in the space J\f corresponds to a 
given kinetic energy K and intersite current J for a par¬ 
ticle on the two-site lattice. This physical parametriza¬ 
tion is universally applicable to lattice systems with any 
number of sites. In the particular two-site case one can 
propose an alternative parameterization the driving field, 
which has an intuitive interpretation in the physical con¬ 
text of the spin-1/2 system. Below we map the space 
M = {K^ J} onto a Bloch sphere and rearrange Eqs. (35) 
accordingly to relate the driving field to a given trajec¬ 
tory in the projective Hilbert space for spin-1/2. 

As a first step we represent the state vector of spin-1/2 
as follows 

\xl)(t)) = e‘^ [cos(6»/2)| t) + e“^sin(6»/2)| ),)] (36) 

where 0{t) and (j){t) are the spherical angles representing 
a point on the Bloch sphere, and /3(t) is an overall phase 
of the wave function. Next, to map the trajectory Q{t) 
to a trajectory on the Bloch we use Eq. (35) and express 


Q in terms of the wave function amplitudes and 

the relative phase (j) 


Q = 20 


l^tP - IV’iP 




(37) 


By substituting from Eq. ( [3^ we relate the com¬ 

plex coordinate Q in the space A/to the spherical coor¬ 
dinates (6>, (j)) on the Bloch sphere 


Q = \ (<^81 


= - 16 sin 0 tan 0 i ^ sin 6> I 


(38) 


This equation gives the required map between the tra¬ 
jectory in the original space of observables to the cor¬ 
responding trajectory of spin-1/2 on the Bl och s phere. 
Finally, by inserting Q of Eq. (38) into Eq. (35c) for B 


(35c) we get a new analytic representation for compo¬ 


(39a) 

(39b) 


nents of the magnetic field 

Bx = (j) tan 0 cos 0 ^ sin 0, 

By = —0 cos 0 + 0 tan 0 sin 0. 

The spherical coordinates (^, 0) determine the wave func¬ 
tion Eq. ( [ 3 ^ up to a common p hase 0 . The phase ^ 
is calculated directly from Eq. ( |35b ) by substituting the 
expressions of \ 2 p^\ and ReQ in terms of 0 and 0, 

1 • 6 

0(t) = 0(0) + - / dt'0tan6> tan-. (40) 

2/0 2 


Equations (36), and (40) solve the problem of re¬ 


constructing the driving field and the wave function from 
a given trajectory on the Bloch sphere. 

Equations (39) demonstrate one subtlety, which 


is very similar to the v-representability problem in 
TD(C)TDFT. Not all trajectories on the Bloch sphere 
are physically reproducible if the driving magnetic field 
is limited to the (x, ^)-plane. For example, it is impos¬ 
sible to drive the system along the equator (7r/2,0(t)) 
with a finite magnetic field because the right hand side 
in Eq. (39) diverges at ^ = 7r/2. Similarly, any trajectory 
which causes a divergence in right hand side of Eq. (39) 


is not S-representable. All physically allowed trajecto¬ 
ries when crossing the equator should approach it in a 
way that 0tan^ stays finite, which translates to the con¬ 
dition 0^0 when 0 ^ tt{2. In other words, a physical 
trajectory, generated on the Bloch sphere by an in-plane 
magnetic field, can cross the line 0 = nj2 only if it is 
perpendicular to that line at the crossing point. This 
is absolutely clear physically because at any instant the 
magnetic filed generates rotation of the spin vector about 
the direction of B. Therefore the initially coplanar to B 
spin is always driven out of the plane. Apparently when 
reconstructing the driving field from a trajectory on the 
Bloch sphere we should take it from a S-representable 
set containing trajectories which either do not touch the 
equator or cross it perpendicularly. 

Now we are ready to present an explicit example of the 
reconstruction. 


Analytically controlled spin flip: design of a quantum NOT 

gate 

To illustrate our inverse engineering formulas we con¬ 
struct a control pulse which does the operation. Initially 
the system is in the ground state corresponding to some 
magnetic field B(0). During the pulse duration r the 
magnetic field is changing and at the end of the pulse 
returns to its initial value B(r) = B(0) while the system 
is driven to the excited state in the field B(0). There¬ 
fore after the pulse the Hamiltonian returns to the initial 
form, but the direction of the spin is reversed. 

For definiteness we assume the initial/final field in the 
x-direction, B(0) = B(r) = xHq. Therefore the initial 
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state is I ^) = l/\/ 2 (| t) + I i))- The target state which 
should be reached at the end of the pulse corresponds to 
another eigenstate of Sx^ that is | ^) = t) ~ I i))- 

On the Bloch sphere the initial and the final (target) 
states correspond to the points (7r/2,0) and ( 7 r/ 2 , 7 r), re¬ 
spectively. 

The first step in constructing the required control pulse 
is find a S-representable trajectory which at t = 0 starts 
at the point (7r/2,0) and arrives to the point ( 7 r/ 2 , 7 r) at 
the time t = r. We note either boundary point belongs 
to the equatorial line. Therefore the trajectory should 
leave the initial point and arrive to the final point along 
the corresponding meridians. To automatically take care 
of the S-representability we introduce a new independent 
variable 7 = (ptanO. The angle 0 in the relevant range 
of 0 < ^ < TT is related to the new variable as follows 


0 = 2 arctan 


\/ (0/7)^ + 1-77 


(41) 


Now the “dangerous” equatorial points correspond to the 
points of the trajectory with = 0 . 

By re-expressing the magnetic field of Eqs. (39) and 


the common phase P of Eq. ( [4Q| ) in terms of 7 and (p we 
find 


Bx = J cos (p — 


dMh) 
{HiY +1 


sin(/), 


By — 


dMh) 

{HiY +1 


COS (/) + 7 sin ( 


and 


=/3(o) + dt' jp +1 -^ 


(42a) 

(42b) 

(43) 


Now we need to find two functions ^(t) and p{t) which 
will do the required job. The first obvious set of condi¬ 
tions for p{t) is 

0(0) = 0, 0(t) = TT, and 0(0) = 0(t) = 0. (44) 


It follows from Eq. (42) that the requirement B(0) = 
B(r) = 5ox will be fulfilled if 

7 ( 0 ) = -i{t) = Bo, (45) 

and dt{(p/"i) = 0 at the boundary points, t = 0,r. The 
latter condition is satisfied if at t = 0, r the second deriva¬ 
tive of 0(t) vanishes 


7 o) = ^{t) = 0 . 


(46) 


In addition we have to make sure that the ratio 0/7 is 
finite for all 0 < t < r. 

As an example we suggest the following ^(t) and 0(t) 
which fulfill all above conditions 


l{t) = ^-Ba;o(l - 2i/r) [(1 - 2t/Tf + 3 ], (47a) 

/X Trt 1 ,47rt, /. X 

ifit) = --- sin(—). (47b) 



FIG. 5. Trajectory of the spin state on the Bloch sphere 
generated by Eqs. ( [47| for r = 12 and Bq — 1. The trajectory 
leaves ground state | and arrives to the excited state | -^) 
perpendicularly (0 = 0) to the equatorial line. 


Here 7 (t) is a smooth monotonically decreasing function 
antisymmetric with respect to the point t = r/ 2 . It 
goes from Bq to —Bq and crosses zero at the middle 
of the pulse. As an extra condition we required that 
7 ( 0 ) = 7 (r) = 0 which allows to smoothly continue the 
driving field be yond the interval 0 < t < r. The func¬ 
tion 0(t) in Eq. (47b) increases monotonically from 0 to tt 
and has zero first and second derivatives at the boundary 
points t = 0 and t = r, and at t = r/2. The derivative 
0(t) = ^sin^(^) is symmetric with respect to the 
middle point t = r/ 2 . 

The corresponding trajectory on the Bloch-sphere is 
shown in Eig. [^for Bq = 1 and r = 12. The trajectory 
starts from the state | on the equator, goes to the up¬ 
per hemisphere, then at t = r /2 it crosses the equatorial 
line at the point ( 7 r/ 2 , 7 r/ 2 ) and reaches the final state 
I -^) from the lower hemisphere. Because of the special 
symmetry of the generating functions the trajectory has 
a central symmetry with respect to the middle point. 

Substituting 7 (t) and p(t) into Eqs. ( [4^ we find mag¬ 
netic field generating this dynamics. In Eig. [^we plot the 
path of the time-dependent magnetic field vector in the 
x^-plane. Each dot represents the magnetic field vector 
at integer times from 0 to 12 . 

The corresponding wave function is given by Eq. (36) 
with the common phase 0(t) defined after Eq. (43). Eor 
our particular example we easily find the overall phase at 
the end of the pulse, 0(t) = 0(0) — 7 r/ 2 .l^ Hence by the 
end of the pulse the system, initially in the ground state, 
is transferred to the excited state, and the wave function 
acquires a common phase shift of 7r/2. 

One can check that if the pulse defined by Eqs. (42) and 
(47) is applied to the initial excited state 10(0)) =0 ^), 
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FIG . 6. The trajectory traced out the magnetic field vector 
(42) in the xy-plane. Starting from the initial value B(0) = x 
the magnetic field follows the trajectory clockwise and comes 
back at t = T = 12. Each arrow represents the magnetic field 
vector at integer times t = 0,l,...12. 


the system is driven to the ground state | ^). Therefore 
our pulse can be viewed as a realization of a quantum 
NOT gate which transforms (up to a common phase shift) 

I i —) into I —and | —into | i —). 


V. CONCLUSION 

In conclusion, we proposed and worked out the strat¬ 
egy of using the TD(C)DFT observable-potential maps 
to inverse engineer analytically solvable time-dependent 
Schrodinger equations and thus to design analytic pulses 
for quantum control problems. We considered a number 
of situations in which the TD(C)DFT maps are known 
analytically. In all those cases the basic observables, such 
as the density or the current, can be used for a convenient 
and physically intuitive parametrization of the driving 
potential and the corresponding wave function. As the 
first pedagogical example we considered the control prob¬ 
lem for the standard textbook Schrodinger equation de¬ 


scribing dynamics of a single quantum particle driven by 
a time-dependent electromagnetic field. In the most gen¬ 
eral case the driving vector potential and the solution to 
the Schrodinger equation can be uniquely reconstructed 
from a given current. If the dynamics is restricted to 
be generated by a scalar potential, the latter can be re¬ 
constructed from the dynamics of the density or, equiv¬ 
alently, from the given (irrotational) velocity field. We 
have demonstrated how from the general reconstruction 
formulas one can recover the known analytic solutions 
of the time-dependent Schrodinger equation for a driven 
harmonic oscillator with a time-dependent frequenc^EEl. 

In the second part of this work we applied our general 
strategy to a less obvious problem of quantum dynamics 
on lattices (discrete spaces). Here we used the analytic 
maps for the one-particle generalized lattice-TDCDFli^ 
where the basic observables are the intersite current and 
the kinetic energy. As a first illustration of this setup 
we considered a manipulation of a one particle state on 
a finite tight-binding chain. We constructed an analytic 
driving field that generates a fast reshaping of the ini¬ 
tial ground state density to the ground state density of 
another potential. In the second example we demon¬ 
strated the engineering of analytically solvable two level 
Schrodinger equations describing, in particular, dynam¬ 
ics of spin-1/2 system driven by a time-dependent mag¬ 
netic field. Here we constructed a cyclic analytic control 
pulse which works as a quantum NOT gate, that is, it 
flips the direction of the spin for two basis stats. 

This work further develops and extends pr omisin g ap¬ 
plications of TDDFT to the quantum control.More 
importantly, it connects the ides of the density-potential 
mapping in TDDFT and TDCDFT to a wider range of 
coherent controP^ and the state preparation problems in 
the quantum computin^^^^^. We hope this connection 
will be beneficial for either field. 
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